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1.  Introduction: 

During  the  contract  period  the  principal  investigator  and  his  associates,  Dr.  I.P.D.  DeSilva 
and  Mr.  Jim  McGrath,  worked  on  several  problems  related  to  turbulence  in  stratified  fluids.  The 
results  of  these  studies  are  outlined  below.  Relevant  publications  originated  from  the  P.I.'s  group 
during  the  contract  period  are  listed  at  the  end  of  this  report. 

The  research  conducted  has  two  components.  In  the  first  part,  a  model  that  deals  with 
turbulent  mixing  across  sharp  density  interfaces  subjected  to  homogeneous  turbulence  was 
developed.  The  analysis  was  essentially  based  on  the  rapid-distortion-theory  approach  of  Hunt 
(1984,  J.  Fluid  Mech.,  138,  161),  which  reckon  on  the  empirical  observation  that  the  dissipation 
rate  of  turbulent  kinetic  energy  (£)  above  an  interface  is  independent  on  the  height  z.  The 
homogeneous  turbulence  was  modeled  using  the  classical  von  Karman  spectra  and  the  interfacial 
motions  were  assu taxi  to  be  governed  by  the  linear  internal-wave  equations.  It  was  further 
assumed  that  the  first  mode  of  internal  waves  governs  the  interfacial  motions,  thus  simplifying  the 
interfacial  boundary  conditions;  the  latter  have  been  a  major  difficulty  in  modeling  fluid  motions 
near  sharp  density  interfaces. 

Two  configurations,  namely,  (t)  a  turbulent  fluid  layer  separated  from  a  non-turbulent 
heavy  fluid  layer  by  a  density  interface  (single-sided  stirring)  and  (»)  a  density  interface 
sandwiched  between  two  layers  of  equal  turbulent  intensities  and  length  scales  (double-sided 
stirring),  were  considered.  The  former  is  applicable  to  situations  of  mixed-layer  deepening  in 
stably  stratified  natural  water  bodies  due  to  convective  stirring  or  surface-wave  breaking  whereas 
the  latter  is  important  in  analyzing  the  interfacial  migrations  in  density  step  structures  -  for 
example,  when  multiple  turbulent  layers  are  separated  by  density  interfaces,  as  in  thermohaline 
staircases  in  the  ocean  (e.g.,  Kelley  1987,  J.  Phys.  Ocean.,  17, 1633).  Detailed  measurements  on 
the  interfacial  structure  for  the  case  (i)  have  been  made  by  Hannoun  &  List  (1988,  J.  Fluid  Mach., 
189,  211),  when  the  interface  is  forced  by  oscillating-grid  induced  shear  free  turbulence. 
Corresponding  measurements  for  the  case  (ii)  were  undertaken  in  the  second  part  of  our  study. 
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The  existence  of  steady  resonant  modes  at  the  interface  was  identified  when  co  <  Ri/2, 
where  (0  is  the  frequency  of  the  waves,  Ri  =  AbL^/u^  js  the  interfacial  Richardson  number,  Ab  is 
the  interfacial  buoyancy  jump  and  is  the  integral  iengthscale  of  turbulence.  Unsteady  analysis 
show  that  the  growing  waves  efficiently  absorb  energy  from  the  turbulent  region  at  or  near 
to  =  Ri/2,  and  the  transfer  of  energy  to  other  modes  takes  place  thereafter.  The  resonant  internal 
waves  intermittently  break  down  owing  to  local  instabilities  and  dissipate  energy  thus  establishing 
a  quasi- stationary  state  at  the  interface  with  a  finite  internal-wave  energy  level.  Such  breaking 
events  cannot  be  described  by  the  linear  theory  and  modeling  based  on  physical  and  mathematical 
reasoning  was  introduced  to  close  the  problem.  On  the  basis  of  the  assumption  (which  is  verified 
by  the  previous  experimental  observations)  that  the  time- traces  of  the  interfacial  displacement  at  a 
given  location  has  steep  gradients,  the  following  spectral  form  was  proposed  for  the  vertical 


velocity: 


,  i  m®)  forcoxoc 
’•'»  !S«0  for  „<•£  • 


where  I(cd)  is  the  spectral  form  predicted  by  the  linear  theory  and  G)c  is  the  frequency  below  which 
the  non-linear  breaking  of  waves  becomes  important  The  cut  off  frequency  o>c  was  evaluated 
using  the  physical  argument  that  the  ratio  of  the  energy  contained  in  the  vertical  velocity  spectra  in 
non-linear  and  linear  regimes  T  is  constant;  this  assumption  led  to  the  condition  that  ^  =  Ri/2d)c  is 
a  constant,  independent  of  the  Richardson  number. 

It  was  found  that  the  experimental  data  on  vertical  velocity  gives  a  good  agreement  with  the 
predictions  for  case  (i),  if  =  0.5  (corresponding  to  T  =  1.25)  is  chosen.  The  calculations  for 
various  quantities  were  carried  out  using  |i.c  =  0.5  for  both  cases  (i)  and  (ii).  In  calculating  the 
displacement  spectra  using  the  above  modeled  spectra,  the  low  frequency  end  near  to  — >0  behaves 
as  to*2,  and  thus  the  existence  of  a  low  frequency  cut  off,  dictated  by  the  lowest  frequency  of  the 

motion  in  the  mixed  layer,  was  suggested.  The  proposed  form  of  the  displacement  spectra  was 

a.  m  ^(CO.OytO2  for  £0  >  COo 
y*(G),°)=  ,  4  , 

1  ¥33  (too.OytOo  for  co  <  too 
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Based  on  previous  work,  the  non-dimensional  coq  was  selected  as  0.18.  Based  on  the 

above  spectral  forms,  with  ^  =  0.5,  ojq  =0.18  and  the  standard  von  Karman  form  of  spectrum  for 

homogeneous  turbulence,  the  Richardson  number  Ri  dependence  of  non-dimensional  r.m.s. 
, — kl/2  Rii/2 

vertical  velocity  (w2)  and  displacement  1C  )  as  well  as  some  other  measurables  were 
calculated;  the  energy  absorbed  to  the  interface  was  calculated  using  the  unsteady  linear  theory. 
The  principal  predictions  are  as  follows: 

Single-Sided  Stirring  Two-Sided  Stirring 


(w^)1'2  =  1.83  Ri‘1/3 

(w^)  =  3.01  Ri1/3 

n\1/2 

/— vt/2 

(C  )  =  3.24  Ri‘5/6 

(c  )  =  5.40  Ri'5/6 

Apw  =  0.28  Ri1/6 

Apw(0)  =  0.14  Ri'1/6 

However,  it  was  pointed  out  that  although  the  wave-number  frequency  spectrum  Xy(k,<o) 
of  homogeneous  turbulence  used  in  the  present  work  yields  the  correct  form  of  Eulerian  frequency 
spectra  in  the  inertial  subrange,  the  predicted  magnitude  is  somewhat  lower  than  previous 
experimental  observations.  Calculations  performed  with  appropriate  modifications  to  Xy  showed 
that  the  predictions  based  on  T  =  1.0  (corresponding  to  ^  =  0.36)  also  give  good  agreement  with 
the  experimental  data.  The  resulting  predictions  are  as  follows  (note  that  the  unsteady-theory 
prediction  for  Apw  remains  unchanged): 


■ 


A  phenomenological  entrainment  model  was  developed  to  calculate  the  entrainment  rates 
across  density  interfaces.  The  model  predictions  were  found  to  be  in  good  agreement  with  the 
entrainment  data  obtained  in  previous  oscillating-grid  experiments  of  Nokes  (1988,  J.  Fluid 
Mech.,  134,  195)  and  Turner  (1968,  J.  Fluid  Mech.,  33, 639).  Turner  (1968)  conjectured  that  the 
entrainment  rates  for  single  and  double-sided  cases  ought  to  be  the  same.  One  of  the  important 
findings  of  the  present  analysis  is  that  these  two  cases  are  different  and  the  different  entrainment 
rates  observed  by  Turner  (1968)  can  be  real;  The  differences  arise  from  the  fact  that  the  profiles  of 
buoyancy  flux  and  the  energy  absorbed  into  the  interface  are  dissimilar.  The  power- law 
dependencies  of  various  quantities  on  the  Richardson  number  Ri,  however,  were  found  to  be 
similar. 

The  present  analysis  is  complete  than  the  previous  theories  on  turbulent  mixing  across 
interfaces.  Hannoun  and  List  (1988)  used  the  analysis  of  Phillips  (1977,  Dynamics  of  the  Upper 
Ocean,  Cambridge  Univ.  Press),  based  on  the  first  mode  of  internal  waves  at  sharp  interfaces,  to 
interpret  their  results  but  the  latter  does  not  take  into  account  the  effects  of  turbulence  on  the 
interface  that  lead  to  the  resonant  breaking  of  waves.  Mory  (1990,  J.  Fluid  Mech.,  223,  193) 
assumed  that  K-H  instability  resulting  from  the  sloshing  motion  of  the  eddies  on  the  interface  leads 
to  mixing.  The  present  analysis  indicates  that  such  a  mechanism  is  unlikely  to  play  an  important 
role  in  mixing,  except  at  small  Ri. 

In  the  second  part  of  the  investigation,  a  laboratory  experimental  study  aimed  at 
understanding  turbulent  mixing  across  shear-free  density  interfaces  was  conducted  by  employing 
state-of-the-art  experimental  techniques  such  as  laser-induced  fluorescence  and  two-point  laser- 
Doppler  velocimetry.  The  experimental  configuration  consisted  of  a  sharp  density  interface 
sandwiched  between  either  a  turbulent  layer  and  a  non-turbulent  layer  (single-sided  stirring)  or  two 
identical  turbulent  layers  (double-sided  stirring).  The  major  focus  was  on  the  measurement  of  the 
magnitude  of  the  r.m.s.  interfacial  velocity  and  r.m.s.  interfacial  displacement  and  their  spectral 
forms  as  well  as  the  fractal  dimension,  thickness  and  the  local  gradient  Richardson  number  of  the 
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interface.  The  mixing  mechanisms  operative  at  different  Richardson  number  ranges  were  also 
studied  by  employing  extensive  quantitative  flow  visualization. 

The  results  show  that  the  mechanism  of  mixing  differs  for  different  Richardson  number  Ri 
ranges.  When  Ri  <  15  for  single-sided  stirring  or  Ri  <  20  for  double-sided  stirring  the  entrainment 
mechanism  was  found  to  be  the  impingement  of  eddies  on  the  interface  and  splashing  of  the 
heavier  fluid  into  the  mixed  layer.  Wave  activities  were  evident  above  these  Ri  values,  but  the 
waves  were  often  found  to  be  disturbed  and  pierced  through  by  the  eddies.  At  such  intermediate 
Ri,  the  entrainment  appears  to  be  dominated  by  the  shearing  of  interfacial-layer  fluid  by  penetrating 
eddies  with  significant  horizontal  velocities.  For  Ri  >  40  or  so,  the  wave  activities  were  prominent 
for  both  single  and  double-sided  stirring  cases  and  the  breaking  of  these  waves  was  the  dominant 
entrainment  mechanism. 

In  both  cases  the  interfacial  layer  was  found  to  be  thin,  with  h/L^  *  0.3,  and  satisfies  Rii  = 
Abh/ufl2  <  k2,  where  h  is  the  interfacial  thickness.  According  to  Carruthers  and  Hunt  (1994,  J. 
Fluid  Mech.,  Submitted),  the  latter  implies  that  the  interface  should  be  dominated  by  the  first  mode 
of  interfacial  waves;  velocity  measurements  corroborated  this  prediction.  On  the  other  hand  when 
Rii  >  Jt2  different  modes  of  internal  waves  can  be  excited  and  a  velocity  maxima  is  possible  within. 
This  has  been  observed  in  the  experiments  of  Perera  et  al.  (1994,  J.  Fluid  Mech.,  In  Press). 
Evidently,  the  difference  in  the  internal-wave  structure  may  cause  different  mixing  rates  across  the 
interface. 

The  normalized  r.m.s.  turbulent  velocity  and  interfacial-layer  isopycnal  displacements 

_  "2  1/2 

showed  a  variation  with  the  bulk  Richardson  number  as  (w2)j  /uh  =  1.7  Ri  -1/3  and  (£  );  /Lh= 

_  —  1/2 

3.2  Ri  *5/6,  for  the  single-sided  stirring  and  (w2)i  /uh  =  3.8  Ri'1/3  and  (£  )4  /Lh  =  4.5  Rr5/6, 
for  the  double-sided  stirring.  These  results  were  found  to  be  in  good  agreement  with  the 
predictions  of  the  above  model.  Typically,  the  gradient  Richardson  number  Rig  within  the 
interface  was  found  to  be  much  smaller  than  Ri,  and  fluctuates  while  frequently  dropping  below 
1/4  and  becoming  negative  at  times.  This  implies  the  possibility  of  intermittent  breaking  of  waves 
within  the  interface,  which  could  be  observed  during  the  flow  visualization  studies.  Apparently, 
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the  dropping  of  Rig  is  concurrent  with  the  increase  of  wave  amplitude  during  resonant  breaking  of 
waves  (see  Phillips  1977). 

Breaking  waves  produces  a  contorted  (or  fragmented)  interface;  the  iso-pycnal  surfaces 
within  the  interface  may  have  approximate  discontinuities,  and  hence  the  Fourier  transformation  of 
spatial  and  temporal  interfacial  displacement  data  is  expected  to  produce  k*2  (wave  number)  and 
to'2  (frequency)  spectra,  respectively;  the  existence  of  such  a  spectral  regime  was  verified 
experimentally.  Note  that  the  model  described  above  employs  this  mathematical  requirement  for 
discontinuous  functions  to  effect  closure. 

The  fractal  dimension  Dk  of  selected  isopycnal  lines  was  determined  as  a  measure  of  the 
fragmentation  of  the  surface.  Physically,  highly  contorted  surfaces  are  expected  to  yield  larger 
fractal  dimensions  and  vice  versa.  The  fractal  dimensions  of  single  and  double-sided  stirring  were 
found  to  be  the  same  for  Ri  <  40,  in  which  parameter  regime  the  interface  is  dominated  by  the 
inertial  effects.  When  wave  breaking  takes  over  interfacial  mixing  at  Ri  >  40,  the  single  and 
double-sided  stirring  showed  different  fractal  dimensions  for  the  same  Ri.  This  shows  that  the 
interfaces  for  these  two  cases  can  have  structural  dissimilarities  in  the  wave-breaking  regime. 
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